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Abstract 

We consider the problem of determining the energy distribution of quantum states that satisfy 
exponential decay of correlation and product states, with respect to a quantum local hamiltonian 
on a spin lattice. For a quantum state on a D-dimensional lattice that has correlation length tr 
and has average energy e with respect to a given local hamiltonian (with n local terms, each of 
which has norm at most 1), we show that the overlap of this state with eigenspace of energy / is 
at most exp{—{{e — f)‘^cr)^^/n'^Da). This bound holds whenever \e — f \ > 2^^/na. Thus, 
on a one dimensional lattice, the tail of the energy distribution decays exponentially with the 
energy. 

For product states, we improve above result to obtain a Gaussian decay in energy, even 
for quantum spin systems without an underlying lattice structure. Given a product state on a 
collection of spins which has average energy e with respect to a local hamiltonian (with n local 
terms and each local term overlapping with at most m other local terms), we show that the 
overlap of this state with eigenspace of energy / is at most exp(—(e — f)'^/nm?). This bound 
holds whenever \e — f \ > rriy/n. 


1 Introduction 

A question of primary interest for local hamiltonian spin systems is to determine the energy dis¬ 
tribution of natural class of states with respect to a given local hamiltonian. The knowledge of 
energy distribution reveals a lot of information about the nature of the state itself. As we shall 
discuss below, a gaussian distribution of energy can be associated to a product state. On the other 
hand, the well known entangled state (also termed as the ‘cat state’) has energy 

distribution peaking at opposite ends of the spectrum of the hamiltonian: |l)(l|i- Moreover, 

the knowledge of energy distribution plays an important role in the study of thermalization of 
quantum systems. 

The aforementioned question has been well studied in classical setting, important examples of 
which are the Chernoff bound [Che52] and the Central limit theorem (which applies to asymptotic 
regime). Chernoff bound can be informally states as follows. Let Xi, X 2 ■ ■ ■ Xn be independent 
and identically distributed random variables taking values in [0,1] and each having average value 
A. Then Pr(| Ai -|- X 2 ■ ■ ■ + Xn — nA\ > e) < where c is a constant that depends on A. 
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One interpretation of this bound (which was the original motivation in [Che52]) is that it 

def 

provides a recipe for distinguishing between two probability distributions P = J2xP{^) l^)(3^l 

def 

Q = I®)(^I with expectation values A and B respectively. Given n independent samples 

xi,X 2 ■ ■ ■ Xn from either of these distributions, the sum is highly likely to be concentrated 

around nA if the underlying distribution is P and around nB if the underlying distribution is Q. 
A more precise formulation of this idea requires characterizing the trace distance between and 
Q®” as n becomes large, and it has been generalized to the quantum setting in [ACMnT'^'07]. 

Another interpretation of the Chernoff bound, which is the focus of present work, lies in the 
setting of ‘classical’ local Hamiltonian systems. Consider a product state p®”’ on n sites, where 
P Sa;P(^) l^)(^l- H be a 1-local Hamiltonian H such that hi = Yl,x^\^){^\ 

def 

acts non-trivially only on the site i and is same for each site. If A = Th:{phi) is the expectation 
value of p with respect to hi, then nA is average energy of p®”’ with respect to the hamiltonian 
H. Let Yi>riA+e be the projector onto eigenstates of H with energy at least nA -\- e. Then the 
Chernoff bound implies that Tr(/3®”n>„^_|_e) < Thus, the energy distribution of p®” is 

highly concentrated around the average energy nA. 

The energy distribution of a product state for quantum lattice system with inhnitely many 
sites was considered in [GV89] (for translationally invariant systems) and in [HMH04a] (for non- 
translationally invariant systems). These results can be regarded as a generalization of the Central 
limit theorem to quantum systems. A non-asymptotic version of Central limit theorem is the Berry- 
Esseen theorem ([Ber41],[Ess45]), which provides an upper bound on trace distance between energy 
distribution of product state and the normal distribution as a function of lattice size. This upper 
bound goes to zero as lattice size approaches infinity, thus recovering the Central limit theorem. 
Eor quantum states with finite correlation length (which includes product states) on finite sized 
lattice, a quantum version of Berry-Esseen theorem was recently shown to hold in [BC15],[BCG15]. 

These results give a strong indication that states satisfying exponential decay of correlation 
behave similar to product states, even when their energy distributions are measured with respect 
to the eigenspectrum of a non-commuting (but local) hamiltonian. The work [KLW15] goes even 
further to show that non-commuting local hamiltonians themselves have energy spectrum that 
resemble that of a 1-local hamiltonian (although, quite curiously, the same work shows that almost 
all eigenvectors of non-commuting local hamiltonians are highly entangled, in contrast with the 
eigenvectors of 1-local hamiltonians). 

Above mentioned results have added to the growing body of research on general properties of 
local hamiltonian systems, such as the Lieb-Robinson bound [LR72], exponential decay of correla¬ 
tion [Has04], the area laws [Has07, ALV12, AKLV13] and local reversibility [KAAV17], to name 
a few. They have also found several applications in the problem of thermalization of many body 
systems. To start with, one of the first steps towards the problem of locality of temperature ^ was 
taken in [HMH04b]. Crucially using the Central limit theorem obtained in [HMH04a], the authors 
characterized a set of conditions under which a given thermal state of a quantum local hamiltonian 
on a lattice would be close to a tensor product of thermal states on local subsystems on the lattice. 

The work [Cral2] considered the problem of thermalization under random hamiltonians, where 
the hamiltonian was generated via a random unitary on a fixed local hamiltonian H. One of the 
main technical problems in this work was the study of the characteristic function Tr(e*^*^), where 

^which is roughly the problem of assigning a temperature to reduced density matrix of thermal state of a local 
hamiltonian, detailed discussion can be found in [KGK"'"13] 
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2 is the maximally mixed state (which is also a product state on the lattice). The techniques 
were inspired from the proof of Central limit theorem in the works [GV89],[HMH04a], where the 
characteristic function Tr(e*'^*/9) of a product state p had been investigated in detail. 

The quantum version of Berry-Esseen theorem [BCG15] was used to show in [BC15] that Gibbs 
state of a local hamiltonian H at sufficiently high temperature (high enough to ensure a clustering 
of correlation) is indistinguishable, over sufficiently large regions of lattice (that scale sub-linearly 
with lattice size), from the microcanonical ensemble of eigenstates of H which have eigenvalues close 
to the average energy of the Gibbs state. This result bears close resemblance to the Eigenstate Ther- 
malization Hypothesis [Sre94, Deu91], which is a stronger conjecture stating that every eigenstate 
of H with eigenvalue close to the average energy of Gibbs state of H is locally indistinguishable 
from this Gibbs state. 

In present work, we provide further details on energy distribution of states satisfying exponential 
decay of correlation and product states. Our main results can be seen as an analogue of the Chernoff 
bound for quantum lattice systems. 

Our first result concerns states that satisfy exponential decay of correlation on a D-dimensional 
lattice. Well known examples of such states include the ground states of gapped local hamiltonians 
[Has04] and Gibbs state above a finite temperature [KGK'*'13]. In fact, it has been shown in 
[FWB'’“15] that for local hamiltonians exhibiting many body localization and having non-degenerate 
energy spectrum, all eigenvectors satisfy exponential decay of correlation. Thus our result provides 
information about structure of eigenvectors of such hamiltonians and may have applications in the 
phenomena of many body localization. 

Fix a L)-dimensional lattice with spins of arbitrary local dimension sitting on each lattice site. 
Gonsider local hamiltonian H on the lattice with n local interaction terms, such that each local term 
has operator norm at most 1 and its support is a hyper-cube of side length 2k, hence containing 
{2k)^ lattice sites (see Section 2 and Figure 2 for a detailed description of H). 


Theorem 1.1 (Informal). Let p be a quantum state with correlation length a and {H)p Tr^pH) 
be the average energy of p. Let n>j (Tl<j ) be the projection onto subspace which is union of 
eigenspaces of H with eigenvalues > f (< f)- 


For a > ./ 20 (-P‘°g W holds that, 




2Dk _(nofa)£ + i_ 

<(D{a)e~-e 


and Tr{pU<iH)p-na) 


2Dfe (na^cT)P + l 

< 0(a)e~ ■ e 


Formal statement of the theorem is given in Theorem 4.2. Thus in one dimensional spin chain 
(with D = 1), our upper bound decays exponentially with energy, rather than as a gaussian. The 
bound becomes weaker with higher dimensions and is depicted in Figure 1. 

Our second result concerns product states over a collection of spins and does not require any 
underlying lattice arrangement of these spins. It does impose, however, a locality constraint on the 
hamiltonian that acts on these spins. Gonsider a hamiltonian H which is a sum of n terms, each 
term being /c-local (that is, it acts non-trivially on at most k spins) and having operator norm at 
most 1. Let m be the maximum number of neighbours of any local term, where two local terms 
are neighbours if there is a spin on which both act non-trivially (See Section 5 and Figure 3 for 
detailed description of iL). We show the following. 
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Theorem 1.2 (Informal). Consider a product state p with average energy {H)p Tr{pH). Fix 

a real number a > ^ ■ Let n>j (n </) he the projection onto subspace which is union of 

eigenspaces of H with eigenvalues > f (< f)- It holds that 


Lr{pIi>^H)+na) < e 5 


and ^ 

Tr(/3n<(//)_„a) < e“ . 

Formal statement of the theorem is given in Theorem 5.3. The energy distribution is depicted 
in Figure 1. The bound is not only independent of any underlying lattice structure, but is also 
independent of the locality k. This is not surprising, since the quantity n that appears in the 
bound is the number of local terms in H, rather than number of spins on which H acts. Following 
corollary is a restatement of above bound, in terms of number of spins (which we call N) and the 
maximum number of local terms that act on any given spin (which we call g). In the following, we 
also assume that each local term is exactly fc-local. 

Corollary 1.3. Consider a product state p with average energy {H)p Tr^pH). Fix a real number 
e > \J(D{g^kN). Let n>j (tl<f) be the projection onto subspace which is union of eigenspaces of 
H with eigenvalues > f (< f)- It holds that 

Tr(pn>(j,)+J < e 

and 

_P_ 

Tv{pU^^H)-na) < e 

Formal statement of the corollary appears as Corollary 5.4. It shows a gaussian decay for tail 
of energy distribution of product states in the scenario where g, k are constants ^ independent of 

N . 


Related recent works 

A recent work [Kuwl6] has obtained a similar concentration result for product states (Lemma 4 
therein). The key idea is to split the hamiltonian H as H = Hi + H 2 + ..., where each Hi , H 2 ... is 
composed of local terms that are non-overlapping. Then from classical Chernoff bound, the product 
state exhibits a Gaussian decay in energy distribution for each of the hamiltonians Hi, H 2 .... Final 
step (which is also the main argument of the paper) is to combine these tails bounds to obtain 
a final bound for energy distribution with respect to the original hamiltonian H. Unfortunately, 
the techniques do not extend to states satisfying exponential decay of correlation. To establish a 
bound for energy distribution with respect to H, one needs the knowledge of bounds for energy 
distribution with respect to each of the ‘classical hamiltonians’ Hi, H 2 ■ ■ ■■ But even for these 

^An interesting class of local hamiltonian system with constant g, k is the family of hamiltonians defined on 
an expander graph, which has recently been a subject of interest with reference to the Quantum PCP conjecture 
[BH13b, AAV13, AE15, FH14, Eldl5, EH15]. 
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Figure 1: The tail bounds according to Theorem 1.1 (right hand side) and Theorem 1.2 (left hand 
side). The x-axis is energy e and y-axis is the weight Tr(/3ne), where He is projector onto eigenspace 
of H with energy e. Shaded region depicts the part of energy distribution with overall weight at 
most The discontinuous part of the curve is where our results provide no information. We have 
ignored 0(1) constants in the figure. 


classical hamiltonians, no bound is known for states that satisfy exponential decay of correlation 
(apart from Theorem 1.1, to the best of author’s knowledge). We have provided further comparision 
of the bound in [Kuwl6] and Theorem 1.2 in Subsection 5.1. 

A concentration result has been noted in [KAAV17] (Section 5 in this reference) for ground states 
of gapped local hamiltonians on finite dimensional quantum lattice systems, which also exhibit 
exponential decay of correlation ([Has04]). In this work, the probability distribution has been shown 
to be concentrated about the median of the distribution with the weight of the distribution above 
energy e decaying as (/ being the median of the distribution, n being the number 

of local terms in the local hamiltonian and a being the correlation length of the ground state). In 
comparison, we show a concentration about the mean of the distribution for all states satisfying 
exponential decay of correlation. While our bounds are weaker than those of [KAAV17] in higher 
dimensions, it may be noted that we have considered a larger class of states that might possess 
weaker properties than the ground states of gapped local hamiltonians. This behaviour appears 
in the context of area laws as well: ground states of gapped local hamiltonians are known to have 
very good scaling of area laws with correlation length [AKLV13]; whereas a recent observation of 
Hastings [Hasl5] suggests that states satisfying exponential decay of correlation may have much 
weaker dependence of area law with correlation length [BH13a]. 

Our technique and organisation 

The idea behind our approach is straightforward, to compute the moment generating function 
Tr(H'^p) of the energy distribution and then use Markov’s inequality to upper bound the desired 
probability. Without loss of generality, we can assume that H = where tc is a label for 

local terms and {hw)p Tr(/9/i^) = 0. Our key technical contribution is the combinatorial lemma 
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(Lemma 3.1) which answers the following question: if we expand as a sum of product of local 
terms, that is = J2wi W 2 -Wr ... hwr, how many terms hu,^hu ,2 ■ ■ ■ ^wr make non-negligible 

(or non-zero) contribution to the moment generating function? We observe the terms that make 
non-negligible contribution possess a common property: there is no hyj^ which is supported ‘far’ 
from all of hwi, Making the notion of ‘far’ precise, we compute the 

number of such terms in Lemma 3.1 and use it to bound the moment generating function. 

The paper is organized as follows. We state basic facts and describe our physical set-up needed 
for Theorem 1.1 in section 2. We prove our combinatorial lemma in Section 3. In Section 4 we 
prove our bounds for states satisfying exponential decay of correlation. In Section 5, we introduce 
the physical set-up required for Theorem 1.2 and provide the proof of the theorem. This proof also 
requires a variant of the combinatorial lemma (Lemma 3.1) which we prove in Appendix A. We 
conclude in Section 6 and address some questions left open by this work. 

2 Physical set-up and basic facts 

In this section, we introduce the physical-set up required for Theorem 1.1. For simplicity of the 
presentation, we shall assume that the spins are arranged on a square lattice, with a local interaction 
term acting between only those spins that are the vertices of a common ‘unit-hypercube’. We 
shall introduce the notion of a ‘dual lattice’ below, to formally and concisely represent these local 
interactions between the spins. It can be observed that more general local interactions on a square 
lattice can be put in this form by sufficient coarse-graining of lattice sites. The physical set-up for 
Theorem 1.2 is relatively simple, and shall be introduced directly in Section 5. 

Consider a H-dimensional real vector space For a vector v € M^, let Vi represent its i-th 
component. For two vectors v^v' G define the ‘1-norm distance’ as 

II /|| def I /| 

It “ II = T* “ '^i\- 

i 

It satisfies the triangle inequality: given G M'^, we have 

||u — u”|| < ||u — U^ll -|- ||u^ — u”||. 

For brevity, we shall refer to 1-norm distance simply as distance. 

For an integer L > 0, define a lattice as the set of all vectors v G that satisfy the 

following: for all i G {1,2...I1}, it holds that Vi is an integer and 0 < u* < L. Two vectors 
v,v' G Td,l are neighbours if ||u — u'|| = 1. Henceforth, the vectors belonging to shall be 

referred to as sites. 

For each site u G £ d,l, we associate a d-dimensional Hilbert space T-L^ and define the full Hilbert 
space as Ji = L^t- Local hamiltonian system is conveniently represented using the notion 

of dual lattice. Let £d,l be the set of vectors w such that for all i G {1, 2... D}, 0 < Wi < L and 
Wi is a half integer (that is, Wi = k + for k an integer). For a fixed w G and an integer k, 

let S{w, k) be the set of all sites v G T_d,l such that: for all z G {1, 2 ... D}, \vi — Wi\ < A: ^ 

A local hamiltonian on is defined as H = l where h^j is a ‘(2A:-|-2)^-local’ term 

that acts non-trivially only on sites in S{w, k) and acts as identity on rest of the sites. The number 
of sites in S{w, k) is at most {2k -|- 2)^, justifying as a ‘(2A: -|- 2)^-local’ interaction. Following 
the physical motivation, we shall refer to vectors in £d,l as interactions. Figure 2 illustrates the 
notions introduced above for the case when D = 2. 
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Figure 2: A lattice £ 2,4 in dimension D = 2 with L = 4. Black dots represent the lattice sites. Blue 
dots represent the vectors of the dual lattice £ 2 , 4 , which we call interactions. For the interaction w, 
the set 5(r/;,0) is {vi,V 2 ,V 3 ,v^}. The distance between sites v and V 3 is ||u — usH = 3. The distance 
between interactions w and tc' is llti; — tc'll = 3. 


Without loss of generality, we assume that the local terms are positive semi-definite ma¬ 
trices and 11 halloo < 1) where ||.||cxd is the operator norm. Given an operator A, support of A 
(called supp(A)) is the set of sites in on which A acts non-trivially. We define the distance 

between two operators A, B to be the minimum distance between their respective supports, that 
is, min.ygg^pp(^^y.y/gg^pp^^^||u H- 

For a quantum state p G Ti, the reduced density matrix of p on a set T of sites is represented 
as pt- We dehne the average energy of p to be Tr(pff), and represent it as {H)p. For every local 

def 

term hw, let {hyj)p = Tii{ps(^yj'^hw). Then we have 

{H) p = 

W 

A state p gB satisfies {C,lQ,a)— decay of correlation if for any two operators A,B such that 
distance between A, B is I > Iq, \i holds that 

\Tr{pA(^B) - Tt{pA)Tt{pB)\ < C||A||||B||e"-. 


Define IIj to be the projector onto eigenspace of H with eigenvalue (energy) equal to /. Let 
n>j (n<y) be the projection onto the subspace which is union of eigenspaces of H with eigenvalues 
greater (less) than /. The following fact follows from Markov’s inequality. 


Fact 2.1. For every t,a > 0 and r even. 


Tr(/9n>(j^)^+„) < 


Tv{p{H - {H)pY) 

{ay 


Proof. We have {H - {H)pY = E/(/ - which gives Tr{p{H - {H)pY) 

{H)pYTv{pIlf). This implies. 


E/(/- 


Tr{p{H-{H)pY)>{aY 

f>{H)p+a 


□ 
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3 A combinatorial lemma 


In this section, we shall prove a combinatorial lemma, which we shall use in Section 4 to prove 
Theorem 1.1. A slight variant of this lemma shall be proved in Appendix A and used in Section 5 
to prove Theorem 1.2. We recall the definition of the set ■Cd,l from Section 2 and let the number 
of interactions in Ld,l be n. It is easily seen that n = {L — 1)^. 

Fix an integer 1. An ordered set {wi,W 2 ■ ■ -Wr} of r interactions in is said to satisfy a 

property P(^) if the following holds: for all Wi, there exists a Wj such that HtCj — WjW < 1. Let the 
number of such ordered sets be N£){n, r, 1). 

Rest of the section is devoted to the proof of following lemma. 

Lemma 3.1. It holds that ND{n,r,l) < ( 4 { 4 l)^nr) 2 . 

We start with the following definition that we shall extensively use. 

Definition 3.2. A selection is an ordered set {(6i,xi), ( 62 , 3 ^ 2 ) ■ • • {br,Xr)}, where 6* € {0,1} and 
Xi S that satisfies the following constraints: 

1. If bi = 0, then Xi can be any interaction in Td,l and if bi = 1, Xi has to satisfy \\xi — Xj\\ < 2-1 
for some j < i. 

2. Number of i for which = 0 is at most |. 

We show the following lemma from which the proof of Lemma 3.1 shall follow immediately. 

Lemma 3.3. Every ordered set {wi,W 2 ■ ■ ■ Wr} that satisfies property P(/) can be mapped to a 
selection in such a way that for any two distinct sets satisfying P{1), the corresponding selections 
are distinct. 

Proof. We assign a selection to an ordered set {vi,V 2 . ■ - Vr} satisfying P(/) using the algorithm 
below. 


Initialization 

• Set i = \ and bi = 0,Xi = Wi. 

• While {i < r), do: 

• z —>■ i + 1. Set Xi = Wi- Set 5^ = 0 if there is no j < z such that ||z«i — Wj\\ < 1. Else set bi = 1. 

• End while. 

Pointer creation 

• Define a relation i? : {1, 2 ... r} —>■ {0,1, 2 ... r} as follows. 

• Set z = 1. While (z < r), do: 

• libi = 1, set R{i) = 0. If bi = 0, find the smallest j > i such that bj = 1 and \\xj — Xi\\ < I (such a j 
exists due to property P(0)- Set R{i) = j. Set z —>■ z + 1. 

• End while. 

Update 

• Let § be the set of all subsets of {1,2,... r} which have cardinality at least 2. 

• For each element S G §, do: 

• Let s be the cardinality of S and Zi, Z 2 ... z® be its elements arranged in increasing order. If it holds 
that bi^ = bi^ = .. .bi^ =0 and R{ii) = i?(z 2 ) = ... R{is) > 0: set bi^ = bi^ = .. .bi^ = 1. 

• End For. 


We show that above algorithm terminates and assigns a selection to each ordered set satisfying 
property P(/). 

1. Consider the running of algorithm during the step Initialization. Condition 1 of a selection 
holds: for every i for which there is a j < z such that \\xi — Xj\\ < I, we have set bi = 1. But 
we haven’t constructed a selection yet, since condition 2 may not be satisfied. 

2. After the step Pointer creation, it may be possible that there exist indices zi,Z 2 ■ - -is (for 
some s < r) such that = bi^ = ■ ■-bi^ = 0, R{ii) = i?(z 2 ) = ■■■R{is) > 0 and is > 
is-i > ...zi. In this case, we find using triangle inequality that Wwi^ — zujj| < Wwi^ — 
WR(i 2 )\\ + \\wR(i^) - zuqll = ||zui 2 - wji{i^)\\ + ||zu/j(i^) - zz;q|| < 21. Similarly, ||zui 3 - zuqH < 
21, |zui4 - zujj <21,... \ wi^ - Will < 2L 

Thus, the step Update sets bi^ = = ... bi^ = 1, recognizing the fact that each of the points 

Wi^, Wi^ ... Wi^ are at a lattice distance of at most 21 from zuj^. This ensures that condition 1 
of selection is still satisfied. 

3. After the step Update terminates, condition 2 of selection is now satisfied as well. We now 

have that for every z with bi = 0, there is no other i! such that R{i) = R{i!^ and bi = = 0. 

Furthermore, = 1- Thus, number of z with = 0 is at most as large as the number of j 
with bj = 1 . 

Lemma follows as two distinct ordered sets satisfying P(f) are not assigned the same selection. 

□ 
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Now we prove Lemma 3.1. 

Proof of Lemma 3.1. For n < r(4/)^, we clearly have r,/) <rf < ((4/)'^nr)5 < (4(4/)'^nr)i. 

So we assume n > r(4/)^. We bound the number of selections, which gives the desired upper 
bound on N£){n,r,l) using Lemma 3.3. 

Consider those selections for which number of i such that bi = 0 is u. For each i with = 0, 
number of possible choices of Xi is n. For each i with = 1 , number of possible choices of Xi is at 
most (4Z)^r (as there are at most (4/)-^ points Xj S Ld,l that satisfy \\xi — XjH < 21 for a given Xi 
^). Hence total number of such selections is at most {^)n^{{Al)^rY~^. Since u < |, total number 
of selections is at most 

r—u ^ 

This proves the lemma. □ 

4 Energy distribution of states that satisfy an exponential decay 
of correlation 

Consider a state p that satisfies (C, Iq, a)— decay of correlation and the hamiltonian H = J2weLo l 
where each term hw is (2k + 2)^-local, that is, it acts non-trivially only on sites in S{w,k). Let 
dw '=^ hw — TT:(phw)i. We prove following bound on r-th moment. 

Lemma 4.1. Given the state p that satisfies {C,lo,a)— decay of correlation, it holds that 
Tr{p{H — {H)Y) < (4(4/o + 8 Dk)^nr )2 + Ce~^a ■ { 4 (^^)^nr ^~^^)2 . 

Proof. Consider, 




Tr(p( 5: gwY) = E Tr(p5 Wl 9 w2 ■ ■ ■ 9wr) ( 1 ) 

W&Id,L WX,W2...Wr 

For every ordered set {wi,W 2 ... Wr} define the quantity D(wi,W 2 ... wfij maxj(minjyj|u;j — 
Wj\). This is the distance of farthest interaction from rest of the interactions in the ordered set. 

For an integer I > 0, define T(l) as the collection^ of all sets {wi,W 2 ...Wr} that satisfy 
D(wi,W 2 ... Wr) = 1. Now, fix a set {wi,W 2 ... Wr} € T{1). Without loss of generality, suppose that 
Wl is an interaction at the distance I from rest of the interactions. The distance between operator 
Qwi and Qwi, for any i 7 ^ 1, is at least I — 2Dk, as the distance from Wi to any site in S(wi, k) is at 
most Dk. Then from (C, a, Iq)— decay of correlation and the relation TT(pgwi) = 0, it holds that 

l—2Dk l—2Dk 

Th{pgwigw 2 ■■■gwr) < T'^{P9 wi) ■ T:r(pgw 2 ■ ■ ■ 9wr) + Ce s— = Ce s—, 

®This is a very crude upper bound and can be found as follows. The number of non-negative integers {ai, a 2 .. .an} 
such that fli < 2? is at most (21)^ . Thus, number of integers { 01,02 ... an} such that |oi| < 21 is at most 
2°{2l)°. 

^to avoid confusion, we call T(l) a ‘collection’ instead of a ‘set’ 


10 





as long as I — 2Dk > Iq. 

Now, the number of sets in the collection 'T{1) is at most Nij(n,r,l) which is upper bounded 
by (4(4/)'^nr) 2 (Lemma 3.1). Thus we have 


^ Tr{pg^^g^^ ... ^ ^ '~^{P9wi9w 2 ■ ■ ■ 9wr) 

Wl,'W2...'Wr I (wi,W2...Wr)£T{l) 

~ 'y ' ^ ' ^^ip9wi9w2 ■ ■ ■ 9wr) 

l<lo+'2Dk {wi,W2...Wr)G'T{l) 

+ E E ^^{p9wi9w2 • • • 9wr) 

l>lo-\-2Dk (wi,W2-..Wr)G‘T{l) 

< Nnin^rJo+ ‘^Dk) + ^ ND{n,T,l)-Ce 

l'>lQ-\-2Dk 

< {A{AIq +'iDk)^nr)^ + Ce~^ ^ (4(4/)'^nr) 2e' 

l>lo+ 2 Dk 

< ( 4 ( 4^0 + 8Dk)^nr)2 + Ce~^ ( 4 nr) 2 e“^ 

i>i 


( 2 ) 


Now, we evaluate 


\ ^ V D I 

^l2 e ^ < 


I 2 e <^dl = a 2 


4^+1 / 


¥e-^ds<a^-t+\—) 


rP I ^ r D , rP 


l>l 


Using this in Equation (2), we obtain 


E "^^{P9wr9w2 ■ ■ ■ 9wr) < (4(4^0 + 8Dfc)^nr)2 + C'e"?E(4(^)^nr^+^)2 


Wi^W2---'Wr 


□ 


Now we proceed to state Theorem 1.1 formally and provide its proof. 

Theorem 4.2. Consider a quantum state p that satisfies {C,lo,a) — decay of correlation and has 
average energy {H)p. 

For > a > ^range exists) it holds that, 


Tr(pn>(^)^+„a) < (1 + Cae^"'°)e s.{4io+8Pk)n Tr(pn<(j^)^_„a) < (1 + C'cre^-'‘)e ss(4io+8Pk)n _ 

For a > holds that, 

2 1 2 1 
. X , 2£)fc . —( 212— r^)D+l ^ ^ ^ 2Dk . —( ^ —rrl^+l 

Tr(pn>(j^)^+„a) < (l+Co-e )e and Tr(pn<(j:^)^_„a) < (l+Ccre )e 

Proof. Using Fact 2.1 and Lemma 4.1 we have, 

,4(4^0+ 8m)^rr A{^)^r^+\r 


Tl{pU>^H}p+na) < (- 
Consider the following two cases. 


na^ 


-)2 + Ce a{- 


na^ 
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> or equivalently ^ 


Then we set r = 2|'(g^^^y7j)^+i], where [.] denotes the ceiling operation (rounding to the 
nearest larger integer) to obtain 


rr. . ^ N A{AlQ + ^Dk)^r,r ^2Dk 4(^)^r'°+\, 

na^ na^ 

, 1 .(4^0 + 8L>A;)^+^ _D_ r mh 

^ T Dna4 )-»>+Ce . a(-), 

< (1 + Co-• 


The last inequality follows from the assumption: - < 1- 


a < 


/8e(4Zo+8i?fc)0+i 1 ^ /8e(4Zo+8i?fc)0 . i ,i 

V nDa - rmd a > y ^ or equivalently 


Dana? 


8 e( 4 Zo+ 8 i?Zc)^ > 1 and 


8e(4Zo+8DZc)(0+i) — 


< 1 . 


We set r = to obtain 

rr. . rr ^ ,4(4^0 + 8T>A;)^r . ^ 2 m 4(^)^r^+\. 

Tr(pn>(H},+na) < - —2 - )-^+Ce ^ a{ ^ -)2 


na^ 


na^ 


Dana^ 


1 T ink- 1 

<1 n I Tia^ 

< (1 + Cae^)e~ 


D\^ 


r) 


Last inequality follows from the assumption: < 1- 

For second part of the theorem, consider the hamiltonian 

H'= Y. 


Define {H')p Tr(pFf') = n— {H)p. Let Ll^j be the projector onto subspace with eigenvalues of H' 
larger than /. Same analysis as above for H' in place of H, along with the relation = n<„_j 
completes the proof. 

□ 


5 Energy distribution of a product state 

In this section, we introduce the physical set-up for Theorem 1.2 and also provide its proof. We 
shall continue using the notations H and h for the hamiltonian and its local term, as this notation 
is restricted only to this section. 

Consider a collection C of spins, such that a d-dimensional Hilbert space "Hf is associated to 
each spin s G C. Let full Hilbert space % be defined as % = For an integer /c > 0, let 
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be the set of all subsets of C of size at most k. For an integer m > 0, let Wk,m be a subset of Sk 
defined as follows (note that Wk,m is also a set of subsets of C) : for each w G VFfc,m the number 
of w' G Wk,m such that \w' H t(;| > 0 is at most m. For each w G Wk,m, let be the set of all 

w' G Wk,m such that Ira H w'\ > 0. Elements of shall be referred to as neighbours of w. The 

set-up has been depicted in Figure 3. 



Figure 3: Collection of spins (blue dots) with local terms (gray triangles). There is no underlying 
lattice structure. Each local term is 3-local, thus k = 3 and each local term has at most 2 neighbours, 
thus m = 2. The set VF 3^2 in above figure is {wi,W 2 ,W 3 , W 4 ^, w^^wq}, where each Wi is the set of spins 
that form vertices of corresponding triangle. Note that for a fixed k,m (here k = 3,m = 2), there 
can be several choices of the set Wk^m and each choice gives a different hamiltonian H. Neighbours 
of W 4 are N('u; 4 ) = {rcsjU^s}. Each spin in the figure is in the support of at most 2 local terms. 
Thus we have g = 2, where g is dehned in Subsection 5.1. 


Let the hamiltonian H be defined as: 


H= Y. 

where acts non-trivially only on spins in w and acts trivially on rest of the spins. Further, we 
assume that ||hi„||oo < 1- 

The definition of Wk,m thus translates to the assumption that: 

1. Each ‘local term’ acts non-trivially on at most k particle, and hence is A:-local. 

2. For each h^, the number of h^i such that the supports of and h^i intersect, is at most m. 

Let p G H he a product state, that is, p = lis^cPs and support of each ps is exactly the spin s. 
Let the reduced density matrix of /? on a subset T C C of spins be denoted in the usual way as pT- 
We bound the moment function Th:{p{H — {H)pY) for an even r to be chosen later and use it 

to prove Theorem 1.2. Define g^ hw — (h^)pL 
We shall prove the following lemma. 

def 

Lemma 5.1. Let n = |114,m| be the number of local terms. Given the product state p = ligi^cPs, 
it holds that 

^(Pi hwT) < (4m^nr)T 
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Proof. Consider, 


Tr(p( i: 9„)’') = i: Tr(/3fl'u)i5ui2 ■ • • 9wr) (3) 

weWk^m Wl,W2...Wr 

Using Tr{pgu)) = Tr(/9^(7^) = 0, we observe that the term Tr(/3g'u,^g'^^ ... gwr) is non-zero only 
if the ordered set {wi,W 2 .. .Wr} satisfies the following property Q: for every Wi, there exists a Wj 
snch that \wi n Wj\ > 0. In other words, there is a wj G Wk^m such that Wi G 'N{wj). Let number 
of ordered sets {wi,W 2 ■ ■ ■ Wr} that satisfy above property be Nk^min, r). This gives us 


■ft(p( E gw) < A^fc,m(n,r)max^j,^2.„^jTr(p5^jC/u,2 .. (4) 

Setting the trivial bound max^j^^2...«)r|3^(w««i5«J2 ■ • ■ gwr)\ < 1; using Lemma 5.2 below, 
the proof follows. □ 

Remark: For the case of translationally invariant systems, where hu, = h for all w, the bound 
max^i,^„ 2 ...^„^Tr(p 5 r^i 5^2 ... g^r) < 1 can be improved to max.ujuw2...wrT^{pgwi9w2 ■ ■■9wr)\ <(1- 
{h)pY. This gives a minor improvement on the statement of Lemma 5.1. 


Now we prove an upper bound on the quantity iVfc „j(n,r) in the following Lemma. Proof is 
very similar to Lemma 3.1 and is deferred to Appendix A for completeness. 

Lemma 5.2. It holds that Nk^rn{n,r) < (4m^nr)2. 

We restate Theorem 1.2 formally and give its proof below. 


Theorem 5.3. Given the product state p 


number a > 



It holds that 


'with average energy {II)p, eonsider a real 


_ na 

T^^{p^>(H)p+na) < e 3^ 


and 

^^iP^<{H)p-na) < e"3^. 

Proof. Lemma 5.1 gives the following upper bound on r-th moment: 


Tv{p{H - {H)p)Y < (dm^nr)^ 

Using Fact 2.1, we have 


Tr(pn>(j:f)^+„„) < 


(4m^nr) 2 
(na)^ 


AfrYr 


r 

2 


Choosing r = 2^^], we obtain for a > 


m ^ TT \ ^ i 

^"(^n>{H>,+na) - (8g^2.„^2 


2 < e 


4em2 ^ 


(5) 
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For second part of the theorem, consider the hamiltonian 


H'= Y. 


Define {H')p Ti{pH') = n — {H)p. Let be the projector onto subspace with eigenvalues of 
H' larger than /. Same analysis as above for H' in place of H gives 


Tr(pnf 


>{H')p+na) — ® 


This completes the proof since = n<„_j. 


□ 


5.1 Restatement of Theorem 5.3 in terms of number of spins 

We introduce a new parameter that captures the number of local terms that act on any given spin. 
Define 

def 1 def 

Qs = Y 9 = max^fit^, 

where Qs is the maximum number of local terms that act non-trivially on spin s. 

Now, we prove Corollary 1.3. Its formal statement is as follows, where we also assume that each 
local term is exactly /c-local. 

Corollary 5.4. Let the hamiltonian H be such that each term hw has locality equal to k. Let 
N \C\ be the number of spins. Given the product state p = lis&cPs with average energy {H)p, 
consider a real number e > y^eg^kN. It holds that 

Tr(pn>(H>^+e) < e 4^^ 

and 

__sl_ 

Tr{pU<{H)p-e) < e 4^^. 

def 

Proof. We set e = na as the energy with respect to H. Then the bound in Theorem 5.3 can be 
restated as: 

Tr(/3n>(/^)^+£) < e 4enm^ . 

Relation between N and n can be computed as follows. To each local term one can associate 
exactly k spins on which hyj acts non-trivially. On the other hand, to each spin s, one can associate 
at most g local terms that contain s in their support. From the first argument, the number of 
associations is exactly k ■ n, whereas from the second argument, the number of associations is at 
most g ■ N. Thus, g ■ N > k ■ n which implies n < Also, m < k ■ g, since each local term is 
supported on k spins, and each of these spins are in the support of at most g other local terms. 
Collectively we obtain nvo? < Ng^k and our bound takes the form: 

Mp^>{H)p+e) < ( 6 ) 

This completes the proof. □ 
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Above upper bound may be compared to Theorem 7 in [Kuwl6]. In this reference, the notion 
of (^'-extensitivity has been introduced (Definition 2, [Kuwl6]), which is analogous to the locality 
parameter g defined above. It is defined as follows: A local hamiltonian H is ^''-extensive if for 
every spin s, we have J2weWkm-s&w W^^wW < g'- Using this, the following theorem has been shown 
in [Kuwl6]: 

Theorem 5.5 (Informal version of Theorem 7, [Kuwl6]). Given a g'-extensive local hamiltonian 
with locality k, it holds that 

Tv{pU>^H),+e) < Q(l)e^P(- ^jy 

where c is a 0(1) constant that depends only on k,g'. 

We observe that Equation (6) achieves a marginally better bound whenever the norm of each 
local term h^ , that is ||/ii„||, is a constant independent of w. In such a case, g' and g are same 
up to the norm of local terms. In case the normalizations of each local term are different, it is not 
clear how g,g' are related to each other. In such a case Equation (6) and Theorem 5.5 may be 
viewed as complementary results. 

6 Conclusion 

We have shown upper bounds on tail of energy distribution of states that satisfy exponential 
decay of correlation and product states, with respect to a local hamiltonian. Main technical tool 
we use is a combinatorial lemma that gives a non-trivial upper bound on the moments of the 
energy distribution. The results may have applications in the study of thermalization of many 
body quantum systems and also for many body localization, as noted in the Introduction. Main 
questions that we leave open are connected to tightness of our bounds, as we discuss below. 

The bounds presented in Theorem 1.2 can only be improved up to constants, since classical 
Chernoff bound also exhibits a Gaussian decay, which is known to be tight. More interesting 
situation occurs with the bounds presented in Theorem 1.1. In one dimensional spin chain, our 
bound decays exponential with the energy. Eor gapped ground states, this is very similar to the 
behaviour noted in [KAAV17] (Section 5) using completely different techniques. This suggests that 
gapped ground states (such as the ground state of Transverse field Ising model, which is exactly 
solvable) are strong candidates for the study of tightness of above results. Our result for higher 
dimensions appears to be much weaker that those obtained in [KAAV17] (Section 5) for gapped 
ground states, and we expect further improvement using better combinatorial arguments. 

An another interesting question is with respect to Matrix product states (with constant bond 
dimension) which are defined on one dimensional spin chain. It is well known that under reason¬ 
able assumptions (see Section 5.1.1, [Orul4]) Matrix product states satisfy exponential decay of 
correlation. Eurthermore, it has already been shown in [OgalO] that given a Matrix product state 
p, if n is large enough and energy e Ri 0(ra), it holds that Tr(pn>^j:^^ _|_£) < It is a strong 

indication that our bound (which applies for all energies e > 0 {y/n)) may be considerably improved 
for this special, but important, class of states. 
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A Proof of Lemma 5.2 

We repeat most of the proof of Lemma 3.1, making changes wherever necessary. 

We start with the following definition. 

Definition A.l. A selection is an ordered set {{bi,xi), (62, X2) ... {br,Xr)}, where bi E {0,1} and 
Xi E H4,mj that satisfies the following constraints: 

1. If bi = 0, then xi can be any element of Wk,m and if bi = 1, Xi has to satisfy |Af(xj)nK(3:j)| > 0 
for some j < i. 

2. Number of i for which = 0 is at most 

We show the following Lemma, from which the proof of Lemma 3.1 follows easily. 

Lemma A. 2. Every ordered set {wi,W 2 ■ ■ .Wr} that satisfies property Q can be mapped to a selec¬ 
tion in such a way that for any two distinct ordered sets satisfying Q, the corresponding selections 
are distinct. 

Proof. We assign a selection to an ordered set {wi,W 2 .. .Wr] satisfying Q using the algorithm 
below. 
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Initialization 


• Set i = 1 and bi = 0, Xi = Wi. 

• While (i < r), do: 

• i —>■ i + 1. Set Xi = Wi- Set 5^ = 0 if there is no j < i such that Wi S J^{wj). Else set bi = 1. 

• End while. 

Pointer creation 

• Define a relation i? : {1, 2 ... r} —)> {0,1, 2 ... r} as follows. 

• Set i = 1. While (z < r), do: 

• li bi = 1, set R{i) = 0. If bi = 0, find the smallest j > i such that bj = 1 and Xi G 'J^ixj) (such a j 
exists due to property Q). Set R{i) = j. Set z —>■ z + 1. 

• End while. 

Update 

• Let § be the set of all subsets of {1, 2,... r} which have cardinality at least 2. 

• For each element S G §, do: 

• Let s be the cardinality of S and Zi, Z 2 .. • Zs be its elements arranged in increasing order. If it holds 
that bi^ =bi^ = .. .bi^ =0 and R{ii) = i?(z 2 ) = ... R{is) > 0: set bi^ = bi^ = .. .bi^ = 1. 

• End For. 


We show that above algorithm terminates and assigns a selection to each ordered set satisfying 
property Q. 

1. Consider the running of algorithm during the step Initialization. Condition 1 of a selection 
holds: for every i for which there is a j < z such that Xi G iN’(xj), we have set bi = 1. But we 
haven’t constructed a selection yet, since condition 2 may not be satisfied. 

2. After the step Pointer creation, it may be possible that there exist indices zi,Z 2 ■ ■ - is (for 

some s < r) such that 6^^ = bi^ = ... bi^ = 0, R{ii) = R{i 2 ) = • • • R{is) > 0 and is > is-i > 
...ii. In this case, we find that Wi^ G G But 'N{wr-^) = ^{{wr^^), 

which implies that |N(zi) n N(z 2 )| > 0. Similarly, |iN’(zi) n N(z 3 )| >0, ... |^(zi) C 3^f(zs)| > 0. 

Thus, the step Update sets bi^ = 6^3 = ... = 1, recognizing the fact that each of the 

points Wi ^, Wi^ ... Wi^ satisfy the property that the neighbourhood of each of them intersects 
with 'N{wij^). This ensures that condition 1 of selection is still satisfied. 

3. After the step Update terminates, condition 2 of selection is now satisfied as well. We now 
have that for every i with bi = 0, there is no other i! such that R{i) = R{i!^ and bi = bi' = 0. 
Furthermore, = 1- Thus, number of i with = 0 is at most as large as the number of j 
with bj = 1 . 

Lemma follows as two distinct ordered tuples satisfying Q are not assigned the same selection. 

□ 


Now we prove Lemma 5.2. 
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Proof of Lemma 5.2. For n < we clearly have <rp < {m?nr )'2 < (Am'^nr)^. 

So we assume n > m^r. We bound the number of selections, which gives the desired upper 
bound on using Lemma 3.3. 

Consider the collection of those selections, for which number of i such that bi = 0 is u. For each i 
with bi = 0, number of possible choices of Xi is n. For each i with bi = 1, number of possible choices 
of Xi is at most (as there are at most number of Xj € Wk,m that satisfy |[N'(xi) H [N’(xj)| > 0 
for a given Xi). Hence total number of such selections is at most Since u < §, total 

number of selections is at most 


E 

u=0 



mfrY “ < ^ [ jn2(m^r)2 < 2^(m^nr)2. 


w=0 


This proves the lemma. 


□ 
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